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q Abstract. Electrostatic spectrometers utilized in high-resolution /3-spectroscopy 

' ^ studies such as in the Karlsruhe Tritium Neutrino (KATRIN) experiment have to 

£>-, operate with a background level of less than 10~ 2 counts per second. This limit can 

*-£ be exceeded by even a small number of 219 ' 220 Rn atoms being emanated into the 

volume and undergoing a-decay there. In this paper we present a detailed model 
of the underlying background-generating processes via electron emission by internal 
conversion, shake-off and relaxation processes in the atomic shells of the 215 < 216 Po 
l/~) daughters. The model yields electron energy spectra up to 400 keV and electron 

r multiplicities of up to 20 which are compared to experimental data. 

CO 



1 INTRODUCTION 2 

1. Introduction 

The Karlsruhe Tritium Neutrino experiment (KATRIN) is a next generation, large-scale 
tritium /3-decay experiment designed to determine the effective electron (anti-)neutrino 
mass m Pe with a sensitivity of 200 meV (90% C.L.) pp. It is currently being assembled 
by an international collaboration at the Karlsruhe Institute of Technology (KIT) in 
Germany. 

KATRIN will investigate the kinematics of tritium /3-decay with unprecedented 
precision in a narrow region close to the /3-decay endpoint Eq ~ 18.6 keV [2]. It is 
only in this narrow region with almost vanishing neutrino momenta that one can gain 
access to m Ve . An essential pre-requisite to obtain the reference sensitivity of 200 meV 
is a low background level of < 10~ 2 counts per second (cps) in the signal region close to 

Eq. 

The KATRIN setup is described in detail in [lj. It consists of a windowless 
gaseous tritium source providing > 10 11 /3-decays per second, a differential and cryogenic 
pumping section to eliminate the injected tritium molecules from the beam line, as 
well as an electrostatic spectrometer acting as high-pass energy filter of unprecedented 
precision, and finally a position sensitive detector to count transmitted electrons. This 
work is focused on background processes in the large spectrometer section. 

In a previous publication [3J we have reported on measurements with the KATRIN 
pre-spectrometer in a test set-up configuration where a-decays of 219 ' 220 Rn atoms in 
the volume of an electrostatic spectrometer of the MAC-E filter typ^f] [U GE |6] were 
identified as significant source of background. These atoms originate mainly from the 
non-e vapor able getter material which is used as a chemical pump to obtain ultra-high 
vacuum (UHV) conditions of p < 10~ 10 mbar [7j, but also from other auxiliary 
equipment within the spectrometer vessel and from the stainless steel vessel hull itself. 
In particular, we could demonstrate that a single radon a-decay can produce up to 
several thousands of detector hits in the energy region-of-interest over an extended time 
period of up to several hours. This background originates from the emission of electrons 
in the energy range from eV up to several hundreds of keV, which is caused by a variety 
of processes related to the emission of the energetic a-particle as well as the subsequent 
reorganization of the atomic shells. Almost all of these electrons are trapped in the 
sensitive volume of the spectrometer due to the known magnetic bottle characteristic 
of a MAC-E filter [8j |9]. Owing to the excellent UHV conditions in this part of the 
setup, electrons remain trapped over very long time periods, so that they can produce 
secondary electrons via ionization of residual gas molecules. A significant fraction of 
these secondaries can reach the detector, resulting in a background rate exceeding the 
KATRIN design limit of 10~ 2 cps. 

In this paper we describe a detailed model of electron emission processes following 
a-decays of the isotopes 219,220 Rn. In a separate publication [10] we in detail validate 
this model experimentally by making use of precise electron trajectory calculations in a 

| Magnetic Adiabatic Collimation with Electrostatic filter 
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MAC-E filter to describe the initial background investigations reported in [HI HH, [12], as 
well as the more in-depth studies performed in [TOj [13]. Furthermore, we make use of 
the model of this work to derive estimates of the background rates and topologies for the 
final KATRIN set-up in [9] , while an active background reduction technique concerning 
trapped electrons is described in [H] . 

This paper is organized as follows: Section [2] will present in detail the processes 
related to electron emission during and after 219 ' 220 Rn a-decay, namely internal 
conversion, inner shell shake-off, relaxation and atomic shell reorganization. The 
implementation of this model into our simulation software will be discussed in section [3j 
Within section [4] we will outline the importance of the physics model implemented in 
this work in the light of our attendant publications [101 [9], ITS] - 

2. Electron emission accompanying radon «-decay 

An essential design feature of high- resolution tritium /3-spectroscopy by a MAC-E filter is 
an excellent UHV in the pressure range p < 1CT 10 mbar, so that background-generating 
ionization processes of /3-decay electrons during the filter process are minimized. In the 
case of the KATRIN spectrometers, this is achieved by non-evaporable getter (NEG) 
strips totalling a length of 3 km in the main spectrometer and 100 m in the pre-spec- 
trometer. As shown in [3], the large surface of the porous NEG strips gives rise to 
emanation of radon atoms associated with the primordial 232 Th, 235 U and 238 U decay 
chains (see figure [I]). Furthermore, the large stainless steel surface of the spectrometer 
vessel (main spectrometer: 650 m 2 , pre-spectrometer: 25 m 2 ) and auxiliary equipment 
attached to it also contributes to radon emanation due to small quantities of radon 
progenitors contained near the surface. 

Due to its long half-life (ti/2( 222 Rn) = 3.82 d [15]), compared to the pump out 
time of radon in the KATRIN spectrometers (£ pr es P ec ~ 25 s, £ m ainspec ~ 360 s), the 
isotope 222 Rn is essentially being pumped out of the spectrometer before it decays. 
Therefore, its background contribution can be neglected. The short-lived isotopes 219 Rn 
(ti/ 2 = 3.96 s) and 220 Rn (t\/2 = 55.6 s), however, will a-decay uniformly over the 
entire spectrometer volume (Vp re spec = 8.5 m 3 , Vmainspec = 1250 m 3 ) to their respective 
daughter nuclei 215 Po and 216 Po (see figure [I]). The important fact for our investigations 
here is that all a-decays are accompanied by the emission of atomic shell electrons from 
the eV up to the multi-keV scale (the a-particle as well as X-ray fluorescence photons 
are of no interest for our background studies, see [9]). If these electrons are emitted 
into the sensitive volume of the spectrometer, they can contribute significantly to the 
background rate via secondary processes. 

There are various processes which can result in the emission of up to more than 
a dozen electrons in a single a-decay. If the a-decay populates an excited level of the 



daughter nucleus, the process of internal conversion, as described in section |2.1[ will 
result in the emission of electrons with energies of up to several hundreds of keV. Also, 
the a-particle itself can interact with electrons of the inner atomic shells, leading to so 
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Figure 1: Top: In the KATRIN spectrometers, non-equilibrium decay chains lead to 
emanation of the two short-lived radon isotopes 219 Rn and 220 Rn. Bottom: radon 
cu-decay processes inside the spectrometer and subsequent electron emission processes 
resulting from shake-off (both isotopes), conversion (mainly 219 Rn) and shell relaxation 
(following conversion and shake-off processes) . 



called shake-off processes, detailed in section 2/2 Both processes produce vacancies in 
the electron shells. These are successively filled by atomic relaxation processes, which 



are the focus of section 2.3 Finally, the shell reorganization process of outer shell 



electrons is described in 12.41 



2.1. Internal conversion 

In an internal conversion (IC) process the excited level of the daughter nucleus, which 
is populated by the a-decay process, interacts electromagnetically with an inner-shell 
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electron, which thus is emitted from the atom. As the IC process is competing with 
radiative processes, it is only dominant for heavy nuclei (P(IC) oc Z 3 ) [161 Eh as * s 
the case for polonium [Z = 84). In addition, the probability of IC decreases for larger 
transition energies, so it is relevant only for low-lying levels. In our specific case, IC 
processes are thus of importance only for 219 Rn — > 215 Po* decays, where significant 
branching ratios lead to the two excited levels (7/2 + , 271.2 keV and 5/2 + , 401.8 keV) 
shown in figure [lj In case of 220 Rn — y 216 Po decays, the even-even nucleon configuration 
of the 216 Po daughter creates a paucity of low-lying excited states, implying that IC 
processes following a-decays of 220 Rn are exceedingly rare processes. 

In an IC process, an inner-shell electron with binding energy E^ is emitted into the 
continuum with a kinetic energy of 

E\dn = E* — E h , (1) 

where E* denotes the excitation energy of the nucleus. For our specific case of 
2i5,2i6p * daughters, conversion electron energies between about 40 keV and 500 keV 
are observed [HI [19] . 

The total conversion probability amounts to about 3.3% when integrating over all 
electron shells in the case of 215 Po*. The probability is largest for K-shell electrons 
(2%) as they are closest to the nucleus. Table [I] lists the dominant electron emission 
probabilities and electron energies for the decay 219 Rn — ?• 215 Po* [18J. Our model allows 
for consecutive IC processes in case the initial de-excitation process does not result in a 
ground state configuration of the polonium daughter. We also include the rare IC process 
of the decay 220 Rn — > 216 Po* [19] , but its contribution is negligible for the investigations 
in [lOj |9]. As mentioned above, the emission of a conversion electron leaves a vacancy 
in the electron shell, leading to subsequent complex atomic shell relaxation processes, 
which are described in sections 12.31 and 12.41 

2.2. Inner shell shake- off processes 

A nuclear a-decay leads to a perturbation of the atomic shells, as the electrons experience 
the passage of the outgoing a-particle through the atomic orbitals, as well as the sudden 
change AZ = Z' — Z of the Coulomb potential of the nucleus (initial state: Z = 86 for 
radon, final state: Z' = 84 for polonium) [20J. The impact of both processes on inner- 
shell (K, L, M) electrons is different than on outer-shell (N or higher) electrons due 
to the largely different orbital velocities. For inner-shell electrons, the orbital period 
is much larger than the orbital passage time of the a-particle (v a /v e ~ 0.1, with v e : 
electron orbital velocity, v a : a-particle velocity). For outer-shell electrons, this ratio is 
reversed (v e /v a ~ 0.1). Accordingly, inner shell electrons will adjust adiabatically to 
the sudden change of nuclear charge. Outer shell electrons, however, remain 'frozen' in 
their parent atom ground state (6p 6 for radon) and will only slowly rearrange to the 



daughter orbitals (6p 4 for polonium), see section 2.4 

For inner shells, electron shake-off (SO) is caused by the direct collision process [201 
22]. In this case, the a-particle, which has already gained 99% of its final kinetic 
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Table 1: The table gives an overview of the relative probabilities Pi (per a-decay) of 
the dominant IC lines and of the corresponding electron energies E^ n for 219 Rn, as 
measured by [IS] . The electron energy is given by eq. Q and can thus be attributed 
to specific excited levels of energy E* via the known values of the binding energy E\> 
(K: 93.1 keV, L: 16.9 keV, M: 4.1 keV, N: 1 keV). Only electron lines with an emission 
probability larger than 0.05% are given. In our model we incorporate the possibility 
of consecutive IC processes within a single a-decay in case that the 401.8 keV level of 
2i5p Q * j g p p U i a ted and de-excites to the 271.2 keV level. 



E kin [keV] Pi [%] shell E* [keV] 



37.5 


0.4 


K 


130.6 


113.7 


0.13 


L 


130.6 


178.13 


1.27 


K 


271.2 


254.29 


0.74 


L 


271.2 


267.08 


0.19 


M 


271.2 


270.24 


0.064 


NP 


271.2 


308.71 


0.233 


K 


401.8 


384.87 


0.102 


L 


401.8 



energy inside the mean radius of the K-shell, can exchange energy with an electron via 
the Coulomb interaction when passing in the vicinity of the corresponding orbital, and, 
consequently, it can kick out the inner shell electron into the continuum. The decay 
energy is shared between the a-particle and the emitted electron. The latter carries 
only a small fraction, usually of the same order of magnitude as the shell binding energy 
£7b- Therefore, in the adiabatic transition, the shake-off process results in a continuous, 
steeply falling energy spectrum. In this work we use the parameterization of Bang and 
Hansteen [23] to determine the emission probability for a SO electron with a certain 
kinetic energy E shakc : 

N(E shakc )=( f h V. (2) 

The SO probabilities for 210 Po have been measured [2U [251 126] and calculated [HJ [27] 
and are used as a good approximation for the 215 > 216 Po isotopes which are considered 
here. The values, which are listed in table |2j underline the well-known fact that the 
ejection probability strongly increases for higher shells. For the M shell, only the total 
emission probability is listed. In our model, we do, however, consider the 5 subshells 
individually, adapting the corresponding emission probabilities. Since there was no 
experimental data available for the individual subshells, we used an equal distribution 
amongst the subshells as an approximation. 
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Table 2: Shake-off probabilities for electrons from specific inner shells in 210 Po, as 
measured by [21J [25] . 



shell probability 



K [2D 1.6 

LI [25] 5.1 

LII [25] 0.6 

LIU [25] 1.5 

M E5] 1.8 



l(T b 

10" 4 
10" 4 
1(T 4 

io- 2 



2.3. Relaxation following conversion and shake- off processes 

An electron, which is emitted via an IC or SO process, will leave a vacancy in an inner 
shell, as shown schematically in figure [2j As a consequence, the electron structure of the 
atom has to rearrange, thereby releasing binding energy. This can be either in the form 
of a fluorescence photon (radiative transition), which is of no concern for this work, or 
in the form of an Auger electron, if the electron filling the vacancy originates from a 
different shell, or a Coster-Kronig electron, in case it is emitted from a sub-shell of the 
same level (non-radiative transition) [28J. In case of a radiative transition, the initial 
vacancy is transferred to a higher atomic shell, while for non-radiative transitions the 
atomic shells are left with two vacancies. The relaxation processes then propagate up 
to the outermost shell. In heavy atoms such as polonium, large cascades are observed 
when inner-shell vacancies are successively filled by non-radiative transitions ("Auger 
explosions"). Consequently, highly charged polonium ions are created, which cannot be 
neutralized when propagating in the spectrometer UHV environment^} 

An electron emitted in a non-radiative transition will receive a distinct kinetic 
energy. In the example of figure [2} the Auger electron energy can, in a first 
approximation, be determined by 

-Ekin = {Eb,Ll — -^b.Ml) — ^b,M4, (3) 

where E^^ are the binding energies of the involved shells i. In case of a radiative 
transition, the photon would have received the energy difference E 1 = -Eb,Li — -^b,Mi- 
The above approximation neglects two effects [29J: 

• A pair of holes in the atomic orbitals retains interaction energy. 

• The relaxation of the atomic orbitals results in a lowering of the final state energy, 
which alters the ionization energies of electron shells containing holes. 

The Auger electron energies, which are applied in our model, are indeed corrected 
for the aforementioned effects, using the intermediate coupling model of [29]. In the 

§ The highly charged Po-ion will be neutralized when hitting the spectrometer vessel. 
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Figure 2: Sketch of a relaxation process. An inner shell vacancy is filled by an electron 
from an outer shell or a neighboring sub-shell, thereby releasing the corresponding 
binding energy difference. This energy can be transferred to another electron, which is 
then ejected from the atom. Depending on the origin of the electron filling the vacancy, 
the emitted electrons are called Auger electrons or Coster-Kronig electrons. 



case of polonium, relaxation electron energies can reach up to about 93 keV, which 
approximately corresponds to the K-shell binding energy. 

Sudden changes of the atomic potentials occur during vacancy cascade development, 
which can lead to the emission of electrons [30]. However, due to their relatively low 
emission probabilities, this effect is only of minor importance and is neglected in our 
model. 

Furthermore, as the number of vacancies in the atomic shells increases, the electron 
binding energies decrease which can lead to the closure of some Coster-Kronig channels, 
reducing the average charge state of the daughter atom [30]. In fact, here we are 
not interested in obtaining an exact multiplicity distribution of emitted electrons 
(predominantly of exceedingly low energies < 100 eV), because the subsequent ionization 
processes of high-energy IC and inner-shell SO electrons in collisions with residual gas 
will produce several hundred or even thousand secondary electrons due to their magnetic 
storage in the spectrometer. 

2.4- Atomic shell reorganization 

In the above described processes the atomic shell of the polonium daughter is left in 
an excited state, and the de-excitation follows a rather complex scheme involving many 
different pathways within relaxation cascades. If the a-decay process leaves the atomic 
shell unperturbed, or if the SO process of the a-particle involves outer shell (N or higher) 
electrons, the relaxation processes are much less complex. The underlying effect is that 
the outer-shell electron wave function cannot adjust adiabatically to the final state due 
to the fact that the outer-shell electron velocity is much smaller than the a-particle 
velocity. In any case the atomic system will relax to the smaller (Z — 2) nuclear charge 
of the daughter nucleus. 

There is a gradual transition of a-decay processes resulting in a highly excited final 
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state to a configuration where the atomic shell is virtually unperturbed, so that the 
initial state with the shell configuration of radon (6p 6 ) has to adjust to the ground 
state shell configuration of polonium (6p 4 ). In our model we treat both processes in an 
identical manner. 

The change in nuclear charge (Z — 2) in the a-decay 86R n — *• 84?° results in a change 
AE = 37.7 keV in the total binding energy of the atomic electrons if the relativistic 
Hartree-Fock-Slater calculations of [31] are used. AE is composed of the sudden 
energy exchange component AE SU( ± and the much slower rearrangement component 
AEr [22]. As the fast inner electrons can adjust adiabatically to the effective nuclear 
charge reduction by rearranging to daughter orbitals, almost all of AE occurs suddenly 
(AE su <i) and has to be supplied by the outgoing a-particle, which results in an equivalent 
retardation. The remaining small fraction AE-& is retained by the atom as temporary 
excitation energy for the much slower shell rearrangement in the outer shells. We employ 
a scenario where the average atomic rearrangement energy AE-& (Qp e — > 6p 4 ) of about 
250 eV [22] is shared statistically by two electrons from the outermost shells. If their 
kinetic energy is larger than the polonium ionization energy for P-shell electrons (1- 
9 eV), they are emitted into the continuum. This results in a flat energy spectrum of 
low-energy "shell reorganization" electrons. 

As the probability for inner shell SO (see table [2]) and IC (see table [TJ is not 
dominating, the above described atomic shell reorganization (SR) in the ground state 
configuration is the most frequent electron emission process accompanying a-decay. If, 
however, the electron shell in the final state is an excited state, caused either by IC or 
by inner shell SO, we calculate the full atomic relaxation, which will be described in 
detail in section |3~T1 

3. Implementation into the simulation software 

To study the event topologies of electrons from the a-decay of 219 ' 220 Rn atoms, and 
to estimate background rates and characteristics due to their subsequent trapping (for 
details, see [lOj [9]), a detailed code for particle trajectory calculations in the complex 
electromagnetic field configuration of the KATRIN spectrometer is required. This 
challenging task is met by the KATRIN simulation package Kassiopeia [32J, which 
allows to track electrons over long periods of time with machine precision. For the 
purpose of this work a Monte Carlo event generator to describe electron emission 



following 219 ' 220 Rn a-decay was developed and is described in section 3.1 Section 3.2 
then gives a short overview of the output of this generator. 

3.1. The radon event generator 

The detailed physical model for signal events and background processes is implemented 
into the Kassiopeia package via MC-based event generators. For the investigations 
of this paper, a radon background generator was developed to describe the processes 
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accompanying the initial radon a-decay, such as internal conversion (IC), shake-off 
(SO), relaxation (RX) and shell reorganization (SR) which are described in the previous 
section |2j 

Figure [3] shows a flowchart of the radon event generator. 
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Figure 3: Event generator flowchart: After initialization of the generator, the different 
physical processes (SO: shake-off, IC: internal conversion, RX: relaxation, SR: shell 
reorganization), represented as solid boxes, are processed according to the model, which 
is described in more detail in the main text. The user has the possibility to configure the 
generator, e.g. turn on or off certain processes to study specific aspects. Corresponding 
decision points are given in dashed boxes. 



The simulation can be configured by the user to study the impact of different 
processes on the background. The following choices are available (options in brackets): 

• activate/deactivate individual physical processes (SO, IC, SR, RX) 

• enforce physical processes (SO, IC) 

• select radon isotope (219, 220) 

During initialization, all data files required for the computation of the various physical 
processes for a specific isotope are read in. Enforcing SO and IC processes can be useful 
because they are rather rare (up to few %). If this option is enabled, it is assured that 
the according processes are executed within every generated event by scaling up the 
emission probabilities of the individual shells until their sum totals 100%. 

The first physical process to be carried out (if activated by the user) is the SO 
process, as it is directly caused by the passage of the a-particle through the atomic 
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shell. At first a random number is generated by the ROOT TRandom3 routine, which 
is based on the Mersenne Twistor algorithm [33]. The SO subroutine then uses the 
generated random number to initiate an SO process with the corresponding probabilities 
for the individual (sub-) shells. Consequently, this can in some rare cases result in the 
emission of multiple SO electrons [20j. For the determination of the SO electron energy 
the acceptance- rejection method [31] is applied to eq. §2§. 

In case of one or more SO electrons being emitted, the RX process will be executed 
(if activated by the user). In our routine we employ the Monte Carlo technique [35] 
to simulate the highly complex pathways of an initial single vacancy, where a large 
number of intermediate electron shell configurations is being involved. In a first step, 
we use the fluorescence yield w, and the Auger yield a, of the shell i which is under 
investigation to determine the transition type. For K- and L-shell vacancies the data 
of [36], and for M- and N-shell vacancies those of [371 EH] is used. If a radiative transition 
is diced, the vacancy is simply transferred to a higher shell, where the new vacancy is 
determined from the available final states according to their relative probabilities. Non- 
radiative transitions up to and including the M-shell result in two vacancies, while 
several vacancies can be created by N-shell vacancy de-excitation due to super-Coster- 
Kronig transitions, i.e. all transitions happen within the N-shell. The described process 
is repeated until all vacancies reach the outer O- or P-shells or until no further de- 
excitations are energetically possible. As we do not take into account small modifications 
of the energies of electron shells due to the actual relaxation process, the de-excitations 
result in discrete energy lines. 

After the RX process was completed, or if, initially, the SO process was deactivated, 
the IC process is performed (if activated by the user). This specific ordering is justified 
by the fact that shell relaxation processes are completed on a much faster time scale 
(10 -15 s) [39] than internal conversion processes (10~ 12 s) [40J. As in the SO subroutine, 
a random number is used to initiate an IC process with the correct probability. Because 
the excited nucleus can decay into an intermediate energy state instead of the ground 
state, this has to be taken into account by allowing consecutive IC processes (so called 
double conversion [22] ). The interrelated energy levels are marked as such in the input 
file, which allows for reliable bookkeeping of the involved states. The IC electron energy 
depends solely on the decaying nuclear state and the binding energy of the emitted 
electron, resulting in discrete IC lines. 

The final process to be carried out is the SR process. At first, the SR subroutine 
checks if any SO or IC processes occurred previously. If this is the case, the SR process is 
skipped because the Po daughter has already relaxed via the above mentioned processes. 
Otherwise, an unperturbed shell is assumed, which results in the excitation of two 
electrons, statistically sharing the shell reorganization energy of AZ?r = 250 eV (shell 
binding energy to be deducted). These electrons are actually only emitted from the atom 
if their energy exceeds the outer shell binding energy of about 1 eV (first ionization) or 
9 eV (second ionization). 

In the final step of the event generation, all electrons generated by a single a-decay 
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are passed to the particle tracking part of the Kassiopeia simulation software. 
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3.2. Generator output 

Figure [4] shows the energy spectra and energy-dependent emission probabilities as 
obtained with our event generator. The discrete IC and relaxation lines, as well as the 
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Figure 4: Event generator: energy spectra of IC, inner-shell SO, relaxation and SR 
electrons for the case of 219 Rn — > 215 Po (top) and 220 Rn — > 216 Po (bottom) a-decay. SR 
electrons, which originate from unperturbed atomic shell relaxation, are distinguished 
from K-, L- and M-shell SO electrons. The electron energy axis is subdivided into 250 
intervals between 1 eV and 500 keV with logarithmically increasing bin size. 



higher-order potential dependence of the SO spectrum can be clearly identified. The flat 
energy spectrum dominating the low-energy part is originating from our model of SR 
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electrons in case of negligible atomic shell perturbation (6p 6 — > 6p 4 ), as the two electrons 
statistically share an energy of 230-250 eV. Due to their identical nuclear charge, the 
inner-shell SO and SR contributions of the two polonium isotopes are assumed to be 
identical. The SO probability is negligible for the inner K-shell. Therefore, the low- 
energy part of the relaxation spectrum results mainly from L-shell (or higher) SO, and 
hence reaches up to about 17 keV (L-shell binding energy). As stated above, the IC 
process is of importance only in the decay 219 Rn —¥ 215 Po*. In this case, there is a high 
probability for vacancies in the inner K-shell, leading to the high-energy part (up to 
about 90 keV) of the relaxation spectrum. 

The total probability for electron emission by a specific process can be obtained by 
integration over the whole energy spectrum. The corresponding results are summarized 
in table [3] 

Table 3: Electron emission probabilities P per decay, depending on the emission process, 
based on the 219 Rn and 220 Rn event generators of this work. P > 1 implies that more 
than one electron can be emitted per decay. 



process P ( 219 Rn) P ( 220 Rn) 

inner-shell SO 2.08 • 10" 2 2.15- 10~ 2 

IC 3.31 ■ 10 2 5.0 ■ 10~ 5 

relaxation 2.29 ■ 10" 1 6.81 ■ 10" 



2 



SR 1.89 1.96 



3. 3. Initial test of the model 

Due to the complex nature of the response of the atomic shells during and after the 
emission of an a-particle, it is of vital importance to compare the present model with 
independent measurements. Generic parameters for comparison are: 
(i) the final charge state of the daughter atom, because it is highly sensitive to a correct 
description of processes such as atomic relaxation, and 

(ii) the electron energy spectrum in the multi-keV range, which can be estimated from 
the number of secondary electrons produced in the electrostatic spectrometer |10j . 
Figure J5J shows the polonium charge state following 219 Rn and 220 Rn decays as obtained 
with the generator of this work, in comparison to the independent measurement reported 
in [26]. There is good agreement between measured and simulated frequencies of 
occurrence of different 216 Po charge states, which underlines the basic validity of our 
model. We attribute the deviations occurring at charge states > 6 to the complex 
nature of atomic relaxation paths, which can only be approximated for the outer shell 
electrons in case of a large initial perturbation of the atomic shell system. 

As can be seen in figure |5j the majority of events consists of two low energetic 
SR electrons. When comparing the experimental data of [26] with our Monte Carlo 
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Figure 5: Charge states of 216 Po (daughter of 220 Rn) and 215 Po (daughter of 219 Rn), 
as obtained with our generator. The simulation is compared to an independent 
measurement of the 216 Po charge state [26J. The values obtained with the 220 Rn 
generator were normalized to the overall rate of non-zero charge states as the 
experimental precision for the detection of neutral daughter atoms was rather limited 
in [26]. For the 219 Rn results, the same normalization constant was used to emphasize 
the difference between the two isotopes. 



generator we note that the detection of neutral daughter states, as outlined in [26], 
was rather challenging. We thus ascribe the discrepancy to experimental difficulties in 
assessing the efficiency in detecting neutral atoms after a-decay. 

The significant difference between the two simulated isotopes in electron 
multiplicities, and correspondingly in the charge distribution of the daughter ion, is 
due to IC processes in the case of 219 Rn — > 215 Po*. As they are emitted from inner 
shells, highly charged final states result from complex relaxation cascades. 

The second important parameter which is of key importance to validate our model is 
the energy spectrum of the emitted electrons in the multi-keV range. In an electrostatic 
spectrometer of the MAC-E filter type, this parameter cannot be measured directly, 
as electrons are trapped over long periods of time [3] |9] [JO] [11] [12] [13]. However, 
an indirect method to assess the energy of stored multi-keV electrons is to make use 
of their subsequent cooling via ionization of residual gas and to count the number of 
produced secondary electrons in a detector. A single radon a-decay can lead to a 
large number of detector hits N^ et (up to 1500 hits corresponding to a single event 
were observed at the KATRIN pre-spectrometer) . There is a good correlation between 
primary electron energy (shown in fig. |4l) and N^et, which is, however, not strictly 
linear due to competing energy losses by synchrotron radiation and due to non-adiabatic 
effects at higher energies. In fig. [6] we display the number of detector hits following single 
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radon a-decays in the KATRIN pre-spectrometer in an experimental configuration where 
ionizing collisions with residual gas (Ar at p — 2 • 10~ 9 mbar) were maximized at the 
expense of synchrotron losses. The measured spectrum is compared to corresponding 
Monte-Carlo simulations with the radon generator of this work. There is good agreement 
between experimental data and MC simulation, taking into account the limited number 
of radon a-decays (127 events) accumulated over a measuring period of about 500 hours. 
The simulation reproduces the main features of the measured distribution: a large 
number of Rn-events with rather few detector counts, caused by the low-energy plateau 
of the shake-off events, and a steep decrease (tail of the shake-off spectrum) towards 
a flat plateau of very few events featuring a large number of detector hits (caused by 
conversion electrons). 
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Figure 6: Comparison of measured and simulated numbers of detector hits produced by 
individual radon a-decay events within the KATRIN pre-spectrometer. An equivalent 
energy scale can be reconstructed when using average energy losses due to scattering 
and synchrotron radiation [10]. The non-linearity of the energy scale results from 
the decreasing scattering cross section for higher energies in combination with linearly 
increasing synchrotron losses. 



A thorough understanding of radon-induced background is crucial for a successful 
neutrino mass determination with the KATRIN experiment. Therefore, further 
detailed background studies comparing measurements and simulations with different 
experimental conditions were carried out in [10]. A combination of the model of 
this work with precise electron trajectory calculations provides the necessary thorough 
understanding in order to reduce the background level below the required limit. 
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4. Conclusions 

In the course of this work we have developed for the first time a comprehensive 
and detailed model of electron emission processes following the a-decays of the two 
radon isotopes 219 Rn and 220 Rn. These investigations were motivated by our earlier 
observations, reported in jH], of periods with significantly enhanced background rates at 
the KATRIN pre-spectrometer measurements. 

The background model incorporates various processes such as internal conversion, 
shake-off and relaxation of the atomic shells during or after the a-emission. The resulting 
electron energies cover a wide range between a few eV up to several hundred keV, and 
involve highly charged polonium daughter ions. Our model successfully reproduces 
experimentally observed polonium charge multiplicities as well as electron energies in 
the multi-keV range, as shown in this work. Further experimental validation of our 
complete physics model was performed in a separate work [10] , where the background 
behavior observed within test measurements at the KATRIN pre-spectrometer could be 
well described. 

It is only by developing and by validating detailed models of background processes 
that the KATRIN experiment can realize its full physics potential in measuring the 
absolute mass scale of neutrinos. 
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